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The hydrodynamic effects on the fate stage kinetics in 
spinodai decomposition of multicomponent fluids are exam- 
ined using a lattice Boltzmann scheme with stochastic fluc- 
tuations in the fluid and at the interface. In two dimensions, 
the three and four component immiscible fluid mixture (with 
a 1024 2 lattice) behaves like an off-criticai binary fluid with 
an estimated domain growth of £°- 4±03 rather than t 1 ^ 3 as 
previously predicted, showing the significant influence of hy- 
drodynamics. In three dimensions (with a 256 3 lattice), we es- 
timate the growth as £°- 96±005 f or both critical and off-critical 
quenching, in agreement with phenomenological theory. 



It is well known that a binary fluid mixture under- 
goes phase separation if rapidly quenched from a high 
temperature phase to a point in the coexistence region. 
Moreover, when the sizes of the domains are much larger 
than the interfacial thickness, there is only one dominant 
length scale in the system It is accepted that the 
late-time dynamics in a binary alloy or glass in which 
the order parameter is conserved, follows a growth law 
of R(t) ~ t 1 / 3 , where R(t) is the average size of the do- 
mains. This growth law is characteristic of the long range 
diffusion of particles between domains and was first pre- 
dicted for off-critical quenches by Lifshitz and Slyozov 
Methods to carry out simulations of phase segre- 
gating systems with hydrodynamic interactions include 
molecular dynamics (MD), direct numerical simulation 
of time-dependent Ginzburg-Landau equations, cell dy- 
namical systems and more recently, lattice gas fl and 
lattice Boltzmann models MM- Lattice Boltzmann sim- 
ulations have recently been successfully used to study 
spinodai decomposition for critical quenches for a binary 
fluid in the presence Q] and absence of porous media 
|^| . The various simulation techniques have tried to ad- 
dress the question of the growth of single phase domains 
and the scaling properties of the correlation or structure 
functions. The underlying aim has been to find certain 
"universality classes" for first order phase transitions in 
which the growth law is independent of the details of the 
interactions, spatial dimension and the number of com- 
ponents. 

The aim of this work is to determine the growth law 
and scaling properties in two and three dimensions for 
three and four component fluids with equal volume frac- 
tions and compare the results to that for an off-critical 
binary quench with unequal volume fractions. We ob- 
tain results in which the hydrodynamic interactions are 



present as well as absent and show that hydrodynamics 
is important. In previous work on the LB method for 
studying phase segregation [4, 5], fluctuations due to cor- 
relations in the particles have been neglected. We include 
the effects of fluctuations on domain growth in multicom- 
ponent systems by incorporating stochastic fluctuations 
in the fluid stress tensor using a scheme recently pro- 
posed by Ladd || and fluctuations in the color gradient 
to perturb the interface. 

The influence of the number of components on domain 
growth and scaling has been recently examined in two 
dimensions for a Potts model with two and three compo- 
nents using Monte Carlo as well as for a three com- 
ponent fluid using Molecular dynamics (MD) Q simula- 
tions. The Potts model does not include hydrodynamics 
and therefore the result of t 1 ^ for the growth law for two 
and three components is not surprising. The conclusion 
of the MD work for a three component fluid shows that 
the growth exponent is also 1/3. Since a thorough MD 
study for fluids requires a large number of particles and 
has to be run for very long times, we have undertaken 
a lattice Boltzmann approach while also incorporating 
fluctuations. Our results demonstrate that for off-critical 
binary quenches and for three and four component flu- 
ids the growth exponent for late times scales as t 03 
in two dimensions. This is to be compared with i 1 / 3 if 
hydrodynamics were not relevant to phase segregation 
and domain growth in two dimensions. Our three di- 
mensional result is the first simulation for an off-critical 
quench to verify the result of Siggia j|] that is based on a 
phenomenological model of droplet coalescence. We also 
determine that dynamical scaling, which is a consequence 
of the existence of one dominant length scale, is valid in 
two and three dimensions. The multicomponent fluids 
show the same scaling behaviour as an off-critical binary 
quench. 

Lattice gas and lattice Boltzmann (LB) methods can 
be used to study hydrodynamic phase segregation using 
parallel computing techniques fllO| ]. They simulate fluid 
properties, phase segregation and the interface dynam- 
ics simultaneously and allow complex boundaries to be 
handled easily. The methods have been described as pro- 
viding the most promising tools to study flow through 
porous media jy] . Unlike methods such as the use of the 
Langevin equation |l2| that are based on a phenomeno- 
logical model of fluid behaviour and are computationally 
intensive, LB methods simulate hydrodynamic phase seg- 
regation in a natural way without the introduction of ad 
hoc relations between the order parameter fluctuations 
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and the fluid dynamics. The lattice Boltzmann method 
is a discrete, in space and time, microscopic kinetic equa- 
tion description for the evolution of the particle distribu- 
tion function of a fluid. Point particles move along the 
links of a lattice (hexagonal in two dimensions), obey 
certain collision rules, and macroscopically mimic the 
Navier-Stokes equations in certain limits. The LB two- 
phase model we used is a modified version of the immis- 
cible fluid model proposed by Grunau et al. |l3| that is 
based on the original model introduced by Gunstensen et 
al. j^J. The multicomponent LB model in this paper is 
a extension of the lattice gas [[l5| and lattice Boltzmann 
p6| models by Gunstensen and Rothman. 
The local order parameter is defined as t) — 

Etc t) ~ Efc=2 *))> Where fit*, t) iS the diS- 

tribution function for the k — th component (fc > 2) 
of the fluid mixture at site x and time t moving along 
the link in the direction i. n is the number of compo- 
nents and N is the number of velocity directions. Also, 
/,(x, i) = Efc=i /z fe ( x ) *) i s t ne distribution function for 
the total fluid, where i = 0, 1, • • •, 6 represent the veloc- 
ity directions at each site of a hexagonal lattice. The 
state i = corresponds to a portion of the fluid at 
rest. The LB equation for /j(x, f) can be written as 
/f(x + ei,t + 1) = /f(x,i) + where k denotes the 
fluid component, and (Of) = (£lf) c + (Of) 33 is the colli- 
sion operator consisting of a term representing the rate 
of change of fj* due to collisions and a term representing 
the color perturbation. The vectors are the velocity 
vectors along the links of a hexagonal lattice. The form 
of (0^) c is chosen to have a single time relaxation with 
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l/ r (/i fe ^ — fi Cq> )> where r is the characteristic 
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relaxation time and eq ^ is the local equilibrium distri- 
bution fl7}| . The surface tension inducing perturbation 
(Of ) p and the recoloring procedure are chosen appropri- 
ately so that Laplaces's law holds for the model [ [l3[ . 
The local color gradient G(x) is defined by G(x) = 

J2iLi e »^( x + e i)i an( i the surface color pertubation that 
is added to segregate and stabilize the interface is 
{n*)P = A | G cos2{6 l - 9 G ), where 9 t is the angle 
of lattice direction i and 9q — &rct&n(G y /G x ) is the an- 
gle of the local color gradient. The surface tension is 
proportional to ~ At p. Here p = ^ k is the local 
particle density. In the original models [jl4 13 1, the re- 
coloring step always makes the local color gradient to lie 
along the direction perpendicular to the interface. To 
mimic the temperature effect on the interface, a noise is 
introduced to perturb the local color gradient direction 
6q, by assuming that the angles are distributed according 
to a Gaussian distribution about Oq- The variance of the 
angle distribution depends on some local temperature. 

A cause for concern about the use of the LB method for 
studying spinodal decomposition has been the lack of sta- 
tistical correlations in the particles. Previous simulations 
show that droplets form and collapse due to an initial ran- 
dom concentration field. It is believed that the spinodal 
decomposition process is strongly correlated with noise 



in the system. It may be argued that the initial fluc- 
tuations may not be equivalent to the natural noise due 
to the temperature of the system. In order to study the 
effects of fluctuations on domain growth, we have incor- 
porated stochastic fluctuations in the fluid stress tensor 
according to a scheme proposed by Ladd The basic 
idea is that, on length scales and time scales interme- 
diate between molecular and hydrodynamic, thermally- 
induced fluctuations can be reduced to random fluctua- 
tions in the fluxes of conserved variables, for instance, the 
stress tensor. It is thus plausible that in an LB simula- 
tion, molecular fluctuations can be modeled realistically 
on intermediate scales, although the microscopic interac- 
tions are different from the real dynamics obtained from, 
for example, an MD simulation. To incorporate this ef- 
fect, a stochastic term, f!(r,t), representing the ther- 
mally induced fluctuations in the stress tensor, is added 
to the time evolution of the density distribution. That is, 
/iCx + ei, t) = / 4 (x, t) + (Of ) c + (Off + //(x, t), where f[ 
is chosen so that its stress moment is nonzero, while con- 
serving mass and momentum. The random stress fluctu- 
ations are uncorrelated in space and time and are sam- 
pled from a Gaussian distribution. The intensity of the 
random stress represents the magnitude of local temper- 
ature. 

We performed critical quenches with < ip(x) > = 0, 
while < ip(x) >^ for off-critical quenches. The largest 
systems simulated were 1024 2 in two dimensions and 256 3 
in three dimensions. Although we have investigated the 
domain growth and scaling properties for a variety of 
lattice sizes and parameters, we report on the domain 
growth and dynamical scaling properties for only one set 
of parameters. The results obtained with smaller lattices 
and different parameters such as surface tension are con- 
sistent with the data presented here. 

The lattice was initialized with a random distribution 
of the different colored fluids. The growth kinetics is 
characterized through the order parameter correlation 
function G(r,t) —< ip(r)ip(0) > — < ip > 2 , averaged 
over shells of radius r. The domain size, R(t), is then de- 
fined as the first zero of G(r, t) and the Fourier transform 
of G{r,t) is the structure factor S(k,t). 

We first discuss the two dimensional results. The ef- 
fects of fluctuations on domain growth is shown in Fig.l 
for a binary fluid after a critical quench. The three cases 
shown are the growth due to (a) thermally induced fluc- 
tuations as described above (+) (b) fluctuations due to 
initial velocity only with average velocity < u >= and 
< u 2 >^ (□) and (c) both thermal and initial velocity 
fluctuations included (O). If neither (a) or (b) is present, 
no patterns are obtained. Since the energy will decay 
due to dissipation in the system, for very long times the 
systems with initial fluctuations or thermal fluctuations 
should behave differently. However, in the current study, 
up to the times simulated, the results are consistent with 
each other and show that the origin of the fluctuations 
has little bearing on domain growth. Thus, the previous 
use [4,5] of an initial fluctuation in velocity only as the 
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driving mechanism for droplet formation would seem jus- 
tified. In all subsequent simulations in this work we have 
used case (c). 

In Fig. 2 is shown the domain growth for a binary 
fluid after off-critical quenchs for (< i)j{x) >= -| (O) 
and < ip{ x ) >— \ (+)• The early time domain growth 
of i 1 / 3 and the later stage growth of t A±m is clearly 
evident regardless of the average order parameter. We 
can interpret this as long range diffusion of fluid parti- 
cles across growing domains at early time giving way to 
domain growth where hydrodynamic or inertial effects 
become important. In order to understand the effects of 
hydrodynamics in the system, the velocity u of the par- 
ticles is set to zero before the collision step. For such a 
system, the convectivc effect and dissipative mechanisms 
will disappear and the only dynamics left is the diffu- 
sive process. In Fig. 3 is shown the domain growth for 
this diffusive motion. The minority phase forms small 
droplets at first which then grow very slowly into larger 
domains. The growth is very slow and is hampered by 
transient effects which slow down the ordering process. A 
transient regime where the growth scales as f' 23 , close to 
t , can be identified and is likely caused by short range 
diffusion along the boundary of the domains, as explained 
by Mullin fll8|| . This short range diffusion crosses over at 
very late times to a t 1 / 3 behaviour characteristic of long 
range diffusion. Though this has been previously pre- 
dicted theoretically and has been detected using ex- 
trapolation methods, we are not aware of any off-critical 
simulations in which it has been directly observed. 

We expect topological factors to play a significant role 
in the ordering process for three and four component fluid 
mixtures. We observe that the domains are rather com- 
pact and can eventually be hindered from growing, sim- 
ilar to the effects of confinement of a fluid mixture in a 
pore geometry jj] . Figure 4 shows the domain growth for 
three and four component mixtures in the presence of hy- 
drodynamics. For the symmetric three (O) and four (x) 
component mixtures an early t 1 ^ 3 regime can be identi- 
fied which crosses over to a clear f 40± - 02 growth for late 
times where the domain morphology is compact. The 
growth for a three component system with concentrations 
(.2, .2, .6) (not shown here) has similar behaviour. Com- 
parison of Figures 2 and 4 shows that, as expected, a 
ternary and four component symmetric fluid mixture be- 
haves like an off-critical binary mixture with an estimated 
growth law of t A± - 03 for late times. The early time be- 
haviour for (O) and (x) shows the t 1 / 3 growth character- 
istic of long range diffusion. Our results clearly indicate 
that hydrodynamics plays a significant role in the late 
time behaviour, contrary to an earlier MD work [|| which 
we suspect was not run long enough. A Langevin model 
simulation of Farrell and Vails |]19[ on a 100 by 100 grid 
for an off-critical binary quench predicts f 35± - 03 . This 
result has been interpreted as showing i 1 / 3 behaviour. 
However, the unpublished work by Wu et al. [^o| with 
this method indicates an exponent very close to results 



reported here. It appears that large lattice sizes and long 
time simulations are necessary to obtain correct growth 
kinetics by this method. The reason for the similar be- 
haviour of three and four component fluids as the off 
critical binary fluid may be understood in terms of the 
number of contact points of different phases. The physics 
appears to be dominated by the two phase line contacts 
rather than the finite number of three or four phase point 
contacts. 

The three dimensional growth is shown in Fig. 5 for 
256 3 for a critical (O) and off-critical (+) quenched bi- 
nary fluid with < ip(x) >= .5. Both cases show similar 
behaviour. Finite size effects have not allowed us to run 
simulation long eough to estimate the exponents for later 
stage growth. However, if we calculate an effective expo- 
nent |2l| , n e ff, by considering the domain sizes at time 
Mt and t, where M is an integer, then the extropolated 
behaviour for large domain is R(t) ~ t°- 96±0 05 (see the 
insert of Fig. 5). The results are in agreement with Sig- 
gia's prediction M. 

Experimental studies has been inconclusive in the es- 
timates of growth exponents. The measurement of an 
off-critical quench of a simple acid-water mixture gave 
a growth exponent between .32 and .35 p3|, whereas a 
growth rate of ~ .5 was observed in a very careful study of 
an off-critical quenched block copolymer [2ljf| . Transient 
effects and the lack of error analysis in the previous stud- 
ies makes direct comparison with experiments difficult. 
Thus we can only stress the need for more experiments 
on off-critical binary quenches and multicomponent flu- 
ids. 

In conclusion, we have used a lattice Boltzmann model 
with fluctuations in the fluid and interface, to study the 
kinetics of domain growth in two and three dimensions 
for an off-critical quenched binary fluid and three and 
four component fluid mixtures. We find that the domain 
growth scales as t 03 in two dimensions, indicating that 
hydrodynamics is relevant to the kinetics of phase sep- 
aration as seen in critical quench. Our three dimension 
results are in agreement with phenomenological theory 
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time diffusive growth, and ~ t A for the later stage 
growth where inertial effects are important. 

Figure 3. Growth for a binary fluid after an off crit- 
ical quench for < ip(x) >— ^ without hydrody- 
namics. The transient behaviour due to interfacial 
diffusion (~ t - 23 ) eventually leads to the expected 
long range behaviour. 

Figure 4. Growth for a symmetric 3-component system 
(O) (1/3,1/3,1/3) and a symmetric 4-component 
mixture (x) (1/4,1/4,1/4,1/4). 

Figure 5. Three dimensional domain growth R(t) on 
a 256 3 lattice for a critically quenched (O) and 
off-critically quenched (< ip(x) >= .5 (+)) binary 
fluid, showing similar growth dynamics. The ex- 
trapolated growth exponent, n e ff versus 1/R(t) as 
seen in the insert, gives ~ t 96 , close to expected 
growth law t 1 . 



Figure Captions 

Figure 1. The effects of fluctuations on domain growth 
for a binary fluid on a 512 2 lattice after a criti- 
cal quench. The three cases represent growth due 
to (a) thermally induced fluctuations via the stress 
tensor (+) (b) fluctuations due to initial velocity 
perturbation only (□) and (c) both thermal and 
initial velocity fluctuations (O). 

Figure 2. Domain growth for a binary fluid after an 
off-critical quench for < ip(x) >= | (O) and 
< ip(x) >= i (+). The straight lines represent 
fits to < ip(x) >= | of ~ i 1 / 3 , showing the early 



